Simple mean-field theory for a zero-temperature Fermi gas at a Feshbach resonance 



in 
o 
o 

(N 

c 



Juha Javanainen, 1 Marijan Kostrun, 1 Matt Mackie, 2 and Andrew Carmichael 1 

'Department of Physics, University of Connecticut, Storrs, CT 06269-3046 
2 Department of Physics, Temple University, Philadelphia, PA 19122 
(Dated: February 2, 2008) 

We present a simple two-channel mean field theory for a zero-temperature two-component Fermi 
gas in the neighborhood of a Feshbach resonance. Our results agree with recent experiments on the 
bare-molecule fraction as a function of magnetic field [Partridge et at, cond- mat/0505353]. Even in 
this strongly-coupled gas of 6 Li, the experimental results depend on the structure of the molecules 
formed in the Feshbach resonance and, therefore, are not universal. 
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Introduction. -Magnetoassociation creates a molecule 
from a pair of colliding atoms when one of the spins flips 
in the presence of a magnetic field tuned near a Feshbach 
resonance [f, 2]. Nevertheless, the ultracold community 
as a whole has only reluctantly acknowledged the role 
of molecules in Feshbach-resonant interactions, clinging 
instead to an all-atom single-channel model solely charac- 
terized by a tunable s-wave scattering length [3] . In this 
Letter we promote a simple molecule-explicit two-channel 
mean-field theory to replace all-atom models. We find an 
excellent agreement with an experiment [4] that should 
present a severe test for theories of all ilk. 



Our approach is built around two species of fermionic 
atoms that magnetoassociate into bosonic "bare" 
molecules. In the experiments [4], a laser spectroscopy 
probe is applied that detects bare molecules only. Below 
resonance, where thermodynamics favors molecules, the 
challenge for single-channel theories is to produce two 
types of molecules out of the same atoms, bare molecules 
seen by the probe and something else. Instead, our two- 
channel model simply states that bare molecules come 
"dressed" with atom pairs. Second, a typical all-atom 
theory rests on the notion of universality, whereby the 
interactions between atoms can be lumped into a scat- 
tering length and microscopic details of the physics do 
not matter. Our approach also mimics universal behav- 
ior on the atom side of the resonance for strong atom- 
molecule coupling. However, the experiments [4] cover a 
wide magnetic-field range on both sides of the resonance 
and, according to our results, certain molecular physics 
details are needed to satisfactorily match theory and the 
experiments. Universality does not hold despite strong 
interactions. Elementary as our mean-field theory is, it 
neatly bypasses two limitations of all-atom theories. 

Model- As before [5, 6], we study two fermionic species 
of atoms (ckj and c k j,) that magnetoassociate into 
bosonic molecules (6k) within a version of the Bose- Fermi 
model [7] used to describe superconductivity beyond the 
usual [8] weak-coupling limit. The corresponding micro- 



scopic Hamiltonian is 
H 
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Here fie k = h 2 k 2 /2m is the kinetic energy for an atom 
with wave vector k, 8 is the detuning controlled by the 
magnetic field (8 > corresponds to an open dissocia- 
tion channel for the molecules), and «k,k' are the atom- 
molecule coupling coefficients for two atoms with mo- 
menta Kk and frk.' . 

The key approximation is to treat the boson operators 
6k and 6 k as classical conjugate variables. Then the hier- 
archy of the equations of motion for fermion correlation 
functions truncates exactly at the level of pair correla- 
tions. As a technical approximation, we only keep the 
zero momentum bosons, i.e., a uniform molecular con- 
densate. We write fermion occupation numbers and pair 
correlations as a function of free-atom energy He = Ke^ as 
P ( e ) = ( c ki c k|> = ( c kT c kT>' C ( e ) = ( c -kic kT ), and, given 
the invariant atom number N arising from the conserved 
quantity = J^kO^k&k + c k | c kt + c k ^ c ki), define the 
molecular amplitude (3 — y/2/Nbo such that f3 2 is the 
fraction of atoms converted into molecules. The result is 
the equations of motion [5] 
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Here S7 = y Nkq_o is the atom-molecule Rabi frequency, 
and /(e) [with /(0) = I] conveys the energy depen- 
dence of atom-molecule coupling. Since the coupling 
coefficient ko.o is inversely proportional to the square 
root of the quantization volume V , the Rabi frequency 
is proportional to the square root of the invariant den- 
sity p = N/V. Finally, he F = H 2 (3t: 2 p) 2/3 /2m is the 
usual Fermi energy for a two-component gas which, so 
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far, arises simply because the sum over wave vectors is 
replaced by an integral over the frequencies. 

The objective is to find the steady state of Eqs. (2)-(4) 
at zero temperature. These equations admit solutions 
of the form C(e;t) = C{e)e'' 2 ^ t , 0(t) = /Je" 2 ^* and 
P(e;t) = P(e), where hji and 2h[i are recognized as the 
chemical potentials for atoms and molecules. In fact, 
there are too many solutions. We therefore impose the 
condition of maximal 'pairing, whereby atoms only come 
in pairs of kj and — kj.. This is technically the usual BCS 
ansatz, but here it is not so much an ansatz as a natural 
consequence of how a pair of atoms couples to a molecule. 
In mathematical terms we have | C (e) | 2 = P(e) — P 2 {e). 
Equations (2) and (4) then give the counterpart of the 
gap equation in BCS theory, 



5-2n = 
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where A(e) = /(e)/3f2/\/2 is the analog of the pairing 
gap. As we [5, 6, 9] and others have done in both bo- 
son [10] and fermion [11] systems, the term oc l/\/e has 
been added into the integrand manually and the corre- 
sponding (possibly infinite) constant on the left-hand side 
has been absorbed into a "renormalization" of the detun- 
ing S. One more equation is needed, and is found from 
the conservation of atom number N that holds for the 
mean- field equations (2)- (4) as well. We have 
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1. (6) 



So far our approach deals with bare atoms and 
molecules that would emerge if the Feshbach-resonance 
interaction were suddenly turned off. In more practi- 
cal terms, the sample breaks up into bare atoms and 
molecules if the magnetic field is switched quickly enough 
far enough away from the Feshbach resonance. The rate 
of change of the detuning that separates quick from slow 
is | <5 1 ~ fi 2 [5, 9], and fast switching may be a tall order 
for the 834 G Feshbach resonance in 6 Li. Here, though, 
alternative methods to access the bare objects have been 
demonstrated. The resonance deals with spin-1/2 atoms, 
the dressed molecules are triplets, and the atoms have a 
net spin equal to zero. A measurement of the magnetic 
moment per atom [12] should [13] therefore only see the 
bare molecules. A laser spectroscopic probe that couples 
only to the triplet molecules likewise allows for a mea- 
surement of the bare-molecule fraction [4] . 

An alternative angle helps with the interpretation of 
the results. Let us initially take a single zero-momentum 
bare molecule. The available state space is spanned by 
the states with either one molecule at zero momentum 
and no atoms, &J|0), or no molecule and pairs of atoms 
k| and — kj, cj^cLkjO), where |0) is the particle vacuum. 



We are at liberty to write the state vector as 

\m) = ^ mbl\o) + -L E c k (*)4 Ajo) . (?) 

Assuming self-consistently that the C coefficients only 
depend on energy, Ck(i) = C(e;t), the Hamiltonian (1) 
gives the time dependent Schrodinger equations for the 
coefficients (5 and C(e) that coincide with Eqs. (2) 
and (4), except that in the counterpart of Eq. (2) we 
have unity in lieu of the factor [1 — 2P(e)]. Were it not 
for this factor that obviously reflects the Pauli exclusion 
principle, our mean-field theory would represent a collec- 
tion of independent dressed molecules. 

A stationary state of the form (7) is a molecule 
"dressed" with atom pairs [5, 6, 9, 13, 14]. With the 
same renormalization that was applied in Eq. (5), the 
dressed molecule is dissociated for 8 > and has one 
bound state for 8 < 0. The latter is the bound state that 
molecular spectroscopy [15, 16] detects in a dilute gas. 
The effect of the renormalization is to put the position 
of the Feshbach resonance in a dilute gas at the detuning 

5 = 0. 

The properties of the dressed molecule depend on 
the coupling function /(e), which may be obtained in 
principle from molecular structure calculations. Setting 
/(e) = 1 corresponds to a contact interaction for the con- 
version between atoms and molecules. In reality, though, 
the interaction has a finite range. Corresponding to a 
nonzero spatial range, there is a range in energy; at 
high enough atom energies the coupling between two bare 
atoms and a bare molecule must vanish. As before [5, 9], 
we use here the model in which we simply cut off the cou- 
pling at an energy KM and write /(e) = 9{M — e) , where 

6 is the Heaviside unit step function. The binding energy 
2h\ii\ of the dressed molecule for a negative detuning — \8\ 
may be solved from the transcendental equation 



3 / m 

\8\ — 2\u\ v/|mIv^ arctan i / - — - 

4 v y \n\ 



(8) 



with u> = il 4 /ep. While and cf depend on the density, 
in the single-molecule parameter u> the density depen- 
dence duly cancels. Given the binding energy, the bare- 
molecule fraction [3 2 in the dressed-molecule state vector 
is easily found. A substantial molecular component, say, 
(3 2 > 1/2 emerges in the contact interaction case M = oc 
for detunings 8 < — (77r 2 /512) u> <~ —0.1 u> and in the case 
M -C u> for detunings 8 < -(3/4) y/MQ. Until further 
notice all results are for the contact interaction, M = oo. 

Main Features.-ln the limit of weak coupling, £1 <C ep, 
the relevant demarkation points for the detuning 6 are 
(Feshbach resonance) and 2cf', 2 because it takes 
two atoms to make a molecule. For S < the sys- 
tem is a condensate of bare molecules, for 5 > 2ep a 
weak-coupling BCS supcrfluid. For < 5 < 2ep we 
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have an atom-molecule mixture in which all atoms from 
the Fermi sea above the energy H5/2 have been con- 
verted into bare molecules with this same energy, so 
that j3 2 = 1 - [5/(2e F )] 3 / 2 . As has been noted be- 
fore [5, 6, 14, 17, 18], the Fermi sea of atoms blocks the 
decay of the molecules even if the molecules are above 
the dissociation threshold. 

In the intermediate case with ~ ep numerical results 
show a continuous transformation from bare molecules 
to dressed molecules, to a nondescript mixture of atoms 
and molecules, and finally to a weak-coupling BCS as 
the detuning is varied from well below the Feshbach reso- 
nance to well above it. Typical 40 K experiments [19] with 
£1 ~ 4€f [5, 6] belong to this category. In the present ap- 
proach the calculated fraction of bare molecules at the 
Feshbach resonance is 6%, and ca. 8% for non-zero tem- 
perature [6, 14], enough to cast ambiguity on the obser- 
vation [19] of fermionic condensation. 

For the 834 G Feshbach resonance in 6 Li, £1 3> ep. 
The coupling is now so strong that dissociation of above- 
threshold molecules is no longer Pauli blocked [6], which 
may bring universal single-channel models into play. 
Since the term universality comes in many shades, we 
frame the issues somewhat along the lines of Refs. [20- 
22], in the order of increasing degree of universality. 
First, there are single-channel theories where atoms and 
their interactions are the building blocks, and explicit 
molecules are absent. The question is, do single-channel 
and two-channel descriptions give the same results? Sec- 
ond, can the relevant atom-atom or atom-molecule in- 
teractions be rolled into a single parameter, such as a 
scattering length, so that the results of cither theory or 
experiments are independent of the microscopic physics 
of the interactions? Third, given the interaction parame- 
ter, will the system reach a unitarity limit that does not 
even depend on this parameter anymore as the interac- 
tion strength increases? 

Zero-temperature strong-coupling BCS theory is a 
single-channel model characterized entirely by the scat- 
tering length, and gives results that are independent of 
the scattering length when the scattering length tends 
to infinity. In this strong form of universality the BCS 
theory agrees with experiments carried out close to the 
Feshbach resonance [23-25] , although it may or may not 
be quantitatively accurate [26, 27]. Moreover, as may 
be deduced easily from the arguments of Ref. [20], in 
the case f2 » ep and uniformly for all detunings 5 > 0, 
our two-channel mean-field theory and the standard BCS 
theory give the same results. The product of Fermi wave 
number and scattering length, kpa, is the dimensionless 
scaling parameter of the BCS theory. The translation 
into our model reads kpa = — 3nQ 2 /8ep6. 

On the other hand, on the side of large negative de- 
tunings and in the limit of strong interaction, the two- 
channel theory scales with the dressed-molecule param- 
eter 8/uj. This parameter cannot be expressed in terms 



of kpa and experimental constants such as atom mass, 
so that BCS and two-channel theories are not equiva- 
lent. As the normalization equation in the BCS theory 
is Eq. (6) sans the explicit /3 2 , an increasing fraction of 
bare molecules signals an increasing discrepancy between 
single- and two-channel mean-field theories. 

Non-universal Experiments. -Our two-channel theory 
covers all parameter regions, and has different scaling 
properties on different sides of the resonance. It is there- 
fore of particular interest that the Rice group has re- 
cently measured [4] the bare-molecule fraction of a Fermi 
gas over a wide range of magnetic fields across the reso- 
nance. Comparisons between theory and experiment are 
complicated by the fact that the experimental density 
varies substantially between data points taken at differ- 
ent magnetic fields. Density information can be recov- 
ered from a representation of the data in terms of the 
parameter kpa [4], but kpa is not a scaling variable in 
our theory and we cannot use it unambiguously as the 
abscissa when plotting the calculated results. For these 
reasons we use the average of the Fermi energies deduced 
from the quoted values of kpa as our parametrization 
of the density. Since the interaction-strength parameter 
scales weakly with density, fl/ep cx p -1 / 6 , we still expect 
a qualitatively valid comparison with experiments. 

The Fermi energy here is given by ep = 7.46 x 27rkHz. 
The detuning in terms of magnetic field B is 5 — Afi(B — 
Bo)/h, where A^t = 2/ib is the difference in magnetic 
moments between a molecule and an atom pair, and Bo 
is the position of the Feshbach resonance. The stan- 
dard relation between atom-atom scattering length and 
magnetic field is a = a^ [1 + A/(B — B )], where at, is 
the background scattering length and A is the magnetic- 
field width of the Feshbach resonance. For a comparison, 
we need to modify the conversion of parameters between 
BCS and our approach to exclude the background scat- 
tering length, so that kp(a — ai,) = — 3irn 2 /8Sep. The re- 
lation between Fermi energy and wave number, of course, 
is hep = fr?k F / (2m). These observations, together with 
the known resonance parameters [16], give the molecular 
parameter Cu — 3.34 x 27rTHz, enormously large com- 
pared with any other frequency scale in the problem. Fi- 
nally, the Rabi frequency is given by Q/ep = (uj/ep) 1 ^ 4 = 
146, which indicates the strong-coupling regime. 

We plot in Fig. 1 the bare-molecule fraction (3 2 as a 
function of the magnetic field B from the Rice data (cir- 
cles) and from our calculations assuming a contact in- 
teraction with M — oo (dashed line). In the experi- 
ments the conversion into bare molecules extends a few 
hundred Gauss below the resonance, while the contact- 
interaction theory predicts a scale <~ 0.1 hui/\Afj,\ ~ 10 T 
and completely misses a qualitative trend in the experi- 
mental data. The remedy is to invoke a finite range for 
atom-molecule conversion. We vary the parameter M un- 
til what is essentially the relative root-mean square error 
between the experiments and the theory is minimized at 
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FIG. 1: Bare molecule fraction in a 6 Li gas as a function of 
magnetic field from experiments [4] (circles), from our theory 
with a contact interaction for atom-molecule conversion, M = 
co, (dashed line), and from the theory with energy cutoff 
M — 254 x 2n kHz corresponding to a range of the interaction 
of about 1500 ao (solid line). The dotted vertical line marks 
the position of the Feshbach resonance. 

M = 254 x 2-7rkHz. The corresponding theory curve is 
plotted as the solid line. In view of the qualitative as- 
sumptions about the Fermi energy in the experiments 
and the simplistic model for the energy dependence of 
atom-molecule coupling, the agreement between theory 
and experiment is now excellent. 

The magnetic-moment data [12] on the bare-molecule 
fraction covers a much smaller dynamical range than 
optical spectroscopy, and only on the molecule side of 
the resonance. The corresponding theory in Ref. [13] 
is a dressed-molecule argument for 6 < that explic- 
itly builds in the background scattering length ai, = 
— 1405 ao Interestingly, when converted into a length 
via M = h/(ml 2 ), our fitted cutoff M gives t = 1500 a . 

Conclusions. -A finite cutoff (M) indicates that a con- 
tact interaction independent of the details of molecular 
physics is not sufficient to reach even a qualitative agree- 
ment with the experiments. The Rice data manifestly 
contradicts universality. Besides, experimental detection 
of a bare-molecule fraction < (3 2 < 1 within a presum- 
ably molecular gas poses a severe challenge to any single- 
channel theory. A single-channel model can be retrofitted 
to include molecules as bound states of atom pairs, which 
would be the counterpart of our dressed molecules. But 
what, then, are the separate bare molecules? 

We have studied a simple two-channel mean-field the- 
ory for a zero-temperature Fermi gas at a Feshbach res- 
onance. The model applies seamlessly across the reso- 
nance; no mixing or matching to join the BCS and BEC 
regimes is needed, or allowed. Our results agree with re- 
cent measurements of the fraction of bare molecules at 
various magnetic fields [4] over a dynamical range of five 



orders of magnitude. Even though the 6 Li Fermi gas is 
strongly coupled, the experimental results are not univer- 
sal but reflect a finite range of atom-molecule coupling. 
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vided both data in a digital form, and valuable advise. 
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